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Intermittency in relation with 1// noise and stochastic differential equations 
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One of the models of intermittency is on-off intermittency, arising due to time-dependent forcing of a bi- 
furcation parameter through a bifurcation point. For on-off intermittency the power spectral density of the 
time-dependent deviation from the invariant subspace in a low frequency region exhibits l/-\/7 power-law 
noise. Here we investigate a mechanism of intermittency, similar to the on-off intermittency, occurring in 
nonlinear dynamical systems with invariant subspace. In contrast to the on-off intermittency, we consider the 
case where the transverse Lyapunov exponent is zero. We show that for such nonlinear dynamical systems 
the power spectral density of the deviation from the invariant subspace can have 1//^ form in a wide range of 
frequencies. That is, such nonlinear systems exhibit 1// noise. The connection with the stochastic differential 
equations generating 1//^ noise is established and analyzed, as well. 



The phrase "1// noise" refers to the well-known 
empirical fact that in many systems at low fre- 
quencies the noise spectrum exhibits an approx- 
imately 1// shape. Generating mechanisms lead- 
ing to 1//'^ noise are still an open question. Here 
we analyze nonlinear dynamical systems with in- 
variant subspace having the transverse Lyapunov 
exponent equal to zero. In particular, we ex- 
plore nonlinear maps having power-law depen- 
dence on the deviation from the invariant sub- 
space. We demonstrate that such maps can gen- 
erate signals exhibiting 1/ f^ noise and intermit- 
tent behavior. In contrast to known mechanism 
of 1// noise involving Pomeau-Manneville type 
maps, coefficients in the maps we consider are not 
static, similarly as in the maps describing on-off 
intermittency. We relate the nonlinear dynamics 
described by proposed maps to 1// noise mod- 
els based on the nonlinear stochastic differential 
equations. 



I. INTRODUCTION 

Intermittency is an apparently random alternation of 
a signal between a quiescent state and bursts of activ- 
ity. In 1949, Batchelor and Townsend used the word 
intermittency to describe their observations of the patch- 
iness of the fluctuating velocity field in a fully turbulent 
fluid.— Many natural systems display intermittent behav- 
ior, for example, turbulent bursts in otherwise laminar 
fluid flows, sunspot activity, and reversals of the geomag- 
netic field. Well known models of intermittency include 
the three types introduced by Pomcau and Manneville,^ 
as well as crisis-induced intermittency."^ A different vari- 
ety of intermittency was first reported when synchronized 
chaos in a coupled chaotic oscillator system undergoes the 
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instability as the coupling constant is changed^i^. This 
intermittency is now known as on-off intermittencyi^^— 
On-off intermittency appears in nonlinear dynamical sys- 
tems with invariant subspaces, where the dynamics re- 
stricted to the invariant subspace is chaotic and the sys- 
tem is close to a threshold of transverse stability of the 
subspace. The main difference of on-off intermittency 
from other types is in the mechanism of the origin: on- 
off intermittency relies on the time-dependent forcing of 
a bifurcation parameter through a bifurcation point; in 
Pomeau-Manneville intermittency and crisis-induced in- 
termittency the parameters are static. 

It is known that the on-off intermittency exhibits char- 
acteristic statisticsi^ (i) the probability density function 
(PDF) of the magnitude of deviation p from the invariant 
subspace obeys the asymptotic power-law, p~^'^'^ , with a 
small positive exponent 77, (ii) the power spectral den- 
sity (PSD) of the time series {p{t)} in a low- frequency 
region exhibits a power-law l/\/7 dependence, and (iii) 
given an appropriately small threshold pth, the PDF of 
the laminar duration t takes an asymptotic form t^^/'^ 
in a certain wide range of t^^ Since on-off intermittency 
generates signals having f~^ PSD with /? = 1/2, the 
question arises whether a mechanism of intermittency, 
similar to the on-off intermittency, can yield signals hav- 
ing other values of the exponent /3 of PSD, in particular 
j3 = 1. The purpose of this paper is to investigate this 
question. 

Signals having the PSD at low frequencies / of the 
form iS'(/) ^ 1//'' with /3 close to 1 are commonly re- 
ferred to as "1// noise", "1// fluctuations", or "flicker 
noise." Power-law distributions of spectra of signals with 
0.5 < /3 < 1.5, as well as scaling behavior in general, are 
ubiquitous in physics and in many other fields, includ- 
ing natural phenomena, human activities, traffics in com- 
puter networks and financial marketSjiSiii Many models 
and theories of 1// noise are not universal because of the 
assumptions specific to the problem under consideration. 
Recently, the nonlinear stochastic differential equations 
(SDEs) generating signals with 1// noise were obtained 
in Refs.llJ (see also recent papersHill) , starting from the 
point process model of 1// noisei^^ Yet another model of 
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1// noise involves a class of maps generating intermit- 
tent signals. It is possible to generate power-laws and 
l//-noise from simple iterative maps by fine-tuning the 
parameters of the system at the edge of chaos-!^^ where 
the sensitivity to initial conditions of the logistic map is a 
lot milder than in the chaotic regime: the Lyapunov ex- 
ponent is zero and the sensitivity to changes in initial con- 
ditions follows a power-law jiS Manneville-^ first showed 
that, tuned exactly, an iterative function can produce in- 
teresting behavior, power-laws and 1// PSD. This mech- 
anism for 1// noise only works for type-II and type-Ill 
Pomeau-Manneville intermittency^S Intermittency as a 
mechanism of 1// noise continues to attract attention;^ 

In this paper we consider a mechanism of intermit- 
tency, similar to the on-ofF intermittency, occurring in 
nonlinear dynamical systems with invariant subspace. 
In contrast to the on-off intermittency, we consider 
the case where the transverse Lyapunov exponent is 
zero. In recent years there is a growing interest in dy- 
namical systems which are characterized by zero Lya- 
punov exponents, namely, which trajectories diverge 
nonexponentiallyj^ Chaos in such dynamical systems 
is called weak chaos. By relating nonlinear dynamics 
with the 1// noise model based on the nonlinear SDEs 
we show that for such nonlinear dynamical systems the 
power spectral density of the deviation from the invariant 
subspace can have \/ f^ form, i.e., 1// noise in a wide 
range of frequencies. Thus a generalization of the on-off 
intermittency yields a new mechanism of 1//*^ noise with 
the asymptotically power-law PDF, originated from the 
commonly known phenomenon of the intermittency and 
weak chaos. 

This paper is organized as follows: in Sec.|lT]we propose 
a model of intermittency with zero transverse Lyapunov 
exponent and in Sec. IIIII we present some examples of 
nonlinear maps exhibiting 1// noise. To obtain analyt- 
ical expressions of the PDF and PSD of the deviation 
from the invariant subspace, in Sec. IIVI we approximate 
discrete maps with SDEs. Section |V] summarizes our 
findings. 



II. MODEL OF INTERMITTENCY WITH ZERO 
TRANSVERSE LYAPUNOV EXPONENT 

We consider two-dimensional maps having a skew 
product structure!^ 



(-n+l 



F{Xn) , V: 



n-l-1 



^V'^ni Un 



(1) 



The function G has the property G(x, 0) = and, thus, 
2/ = is the invariant subspace, while y is the deviation 
form the invariant subspace. We assume that the dy- 
namics Xn+i = F{xn) in ([IJ restricted to the invariant 
subspace is chaotic. If the transverse Lyapunov exponent 



1 ^"^ 
lim — > In 

n=0 



5G(x„,0) 



dy 



(2) 



along an orbit on the invariant subspace converges and is 
less than zero, then the invariant subspace is transversely 
stable with respect to this orbit. 

In this article we consider the case when dG{x, 0)/dy — 
1 and, consequently, the transverse Lyapunov exponent 
is zero. Furthermore, we will assume that the two terms 
with the lowest powers in the expansion of the function 
G(x, y) in the power series of y have the form 



Gix,y) ^y + g{x)y'' 



(3) 



with T] > 1. This form satisfies the condition 
dG{x,0)/dy — 1. Particularly 77 = 2, however, gener- 
ally J] may be fractional, as well. 

We will consider the case where the function g{x) in 
Eq. ^ is not constant and can acquire both positive and 
negative values. Thus the expansion ^ leads to the the 
map for small values of y„ 



Vn+l = Vri + Zny"^ , 



77 > 1 



(4) 



where z„ = gixn)- It should be noted that when rj — \, 
the map ^ becomes a multiplicative map yn+i = 2/n(l + 
z„), which is essentially the same as the map considered 
in Ref. |7| for modeling of on-ofF intermittency. The map 
(m is similar to Pomeau-Manneville map 



yn+i^yn+ayl (modi) 



(5) 



on the unit interval with one marginally unstable fixed 
point located at y = 0.^ The main difference from the 
map ((H) is that in the Pomeau-Manneville map ([S]) the 
coefficient in the second term is static. 

Let us consider the situation when j/„ > 0. If z„ < 
then the map (j4|) leads to the decrease of the deviation 
from the invariant subspace y = 0, whereas for z„ > the 
deviation y grows. In contrast to systems with nonzero 
transverse Lyapunov exponent, the growth or decrease 
of the deviation is not exponential. In fact, if the second 
term on the right-hand side of Eq. ([4]) is much smaller 
than the first and, consequently, Eq. (j3]) can be approxi- 
mately replaced by the differential equation dy/dt — y^^z, 
the growth or decrease of the deviation y can be de- 
scribed by a g-exponential function with q = rj. The 
g-exponential function, used in the framework of nonex- 
tensive statistical mechanics, ^'^"^^ is defined as 



exp^{x) = [1 + (1 -q)x]\ 



(6) 



where [x]^ = max{a;,0}. Thus, although the Lyapunov 
exponent is zero, the map can be characterized by a 
nonzero g-generalized Lyapunov coefficient.— i^ 

If the average of the variable z is positive, (z) > 0, and 
there is a global mechanism of reinjection, the map @ 
leads to the intermittent behavior. As in on-off intermit- 
tency, the intermittent behavior appears due to the time- 
dependent forcing of a bifurcation parameter through a 
bifurcation point z = 0, thus the behavior described by 
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map ([4]) can be considered as a kind of on-ofF intermit- 
tency. However, on-off intermittency is usually investi- 
gated in dynamical systems with nonzero transverse Lya- 
punov exponent. 

For small durations of the laminar phase, one can ap- 
proximate the map Q replacing ?/„ in the second term 
on the right hand side with initial value j/o- In this case 
Eq. (m describes a random walk with drift. Since the 
average displacement due to the diffusion grows as \Jt 
and the displacement due to drift term is proportional 
to i, for small enough durations t the diffusion is more 
important than the drift. It is known that for the un- 
biased random walk the distribution of the first return 
times has the power-law exponent — 3/2i^ Therefore, for 
small enough durations t one can expect to observe the 
power-law form, t^'^/^, of the PDF of the laminar phase 
durations, the same as in on-off intermittency. 

The first two terms in the expansion ([3]) do not allow to 
determine uniquely the PDF of the deviation y. In order 
to determine PDF of y and PSD of the series {?/«}, we 
need to take into account more terms in the expansion of 
the function G(x, y) in the power series of y. One of the 
possibilities that we will consider is for the third term in 
the expansion to be equal to 72/^''"^ (note, that 2r\—\ > rj 
when ?7 > 1), leading to the map 



Vn+l = Vn + Znyll + 7y, 



2j)-1 



(7) 



Particularly, for 77 = 2, 2?7 - 1 = 3 and Eqs. ® and © 
display simple Taylor expansions. Note, that a mecha- 
nism of reinjection operates at large values of y and does 
not change Eq. Q , written for small values of y close to 
the inveriant subspace. 



A. g-exponential transformation of random walk 

Another example of the function G{x, y) having the 
expansion in the power series of y as in Eq. ^ can be ob- 
tained according to the following consideration: In Ref. [j 
a map of the form 

2/„+l = Wnyn (8) 

was considered as a model of on-off intermittency. In the 
log domain this map transforms to 



^n+l — ^n ~r Zn 



(9) 



where s„ = lnj/„ and z„ = lni(;„. The critical condition 
for the onset of on-off intermittency is the condition for 
unbiased random walk, (z) — 0. One of the reasons for in- 
termittent behavior is highly non-linear relation y„ — e**" 
between s„ and y„. We can expect intermittent behav- 
ior also using other nonlinear functions instead of the 
exponential function. One of the generalizations of the 
exponential function, which corresponds to the differen- 
tial equation dy/ds = y, is the g-exponential function 
^ obeying the equation dy/ds = j/''. Thus, instead of 



Vn = e''" we will consider a relation j/„ — exp (s„), lead- 
ing to a map of the form 

y„+i = exp^(ln^(y„)-|-z„) = (y^"'' + (I-77)z„)T^ , (10) 

where the g-logarithm, defined as2i 



Ing(x) 



.1-9 „ 1 
1-? '■ 



(11) 



is a function inverse to g-exponcntial function. Expand- 
ing the map PH)) in power series of y we get Eq. Q . 

The g-exponential function exp (s„) tends to infinity 
as s„ approaches 1/(77—1) and the variable y„ = exp (s„) 
can be introduced only when s„ does not reach 1/(77— 1). 
This can be achieved by modifying the map ^ for the 
values of s„ close to 1/(77 ~ 1) in order to avoid reaching 
this value. The modification of the map (jH]) changes also 
the map pO)) for large values of ?/„ , not allowing for the 
value of the expression y^^^ -|- (1 — ri)zn to become zero. 



III. NUMERICAL EXAMPLES 

In this Section we present some examples of the map 
^ with the function G{x, y) whose behavior for small 
values of y is described by Eq. (O or Eq. (fTO|) . Let us 
consider the map ([7|) with 77 = 2, 7 = 0.5 when the vari- 
able Zn has the average (z) = 5 x 1G~^ and the variance 
{{z — {z)Y) = 1. The parameters of the map are chosen 
taking into account equations from Sec. IIVI The chosen 
value of the average {z) is close to the critical value for 
the onset of intermittency (z) = and is much smaller 
than the standard deviation of the variable 2„. As a 
mechanism of reinjection we use a reflection at y = 0.5, 
leading to the map 



y„+i = 0.5 - |y„ + z„y^ + 0.5y^ - 0.5| 



(12) 



As a map Xn+i — F(xr, 
driving by a tent map 



in Eq. ([T]) we take the chaotic 




< a;„ < i 
h<x„< 1. 



(13) 



The variable z„ with given average (z) and variance ((z- 
(z))^) can be obtained from Xn using the equation 



{{x-{x)Y) 



(a:„ - (x» + (z) 



(14) 



For the tent map ([T3| the average and the variance are 
(x) = 0.5 and {{x — (x))^) = 1/12, respectively. 

Another example is when the variable z„ acquires only 
two values ±C, with the probabilities p+ and p_, p+ + 
P- — \. In this case the average and the variance of z„ 
are given by the equations(z) — C,{p+ ~ P-) and ((z — 
(z))^) = 4C^p+p_. Expressing the probabilities we get 



_ 1 (z) 

^^ " 2 ^ 2V((z-(z))2)-H(^ 



(15) 
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and 



c = ^{{z-{z)r) + {zy 



(16) 



In particular, if ((z— (z))^) = 1, (z) = 5x 10 ^ thenp+ w 
0.500025, p- w 0.499975, C ~ 1.00000000125. Such two- 
valued variable z„ can be implemented by the following 
map: 



yn+i 



Vn - vie + 0.5y3 , 

0.5 - |y„ + y^C + 0.5y3 



^ a;„ ^p_, 

0.5|, p_ <a;„ ^ 1. 

(17) 

Note, that also in the map (IT7|) we use a reflection at 

y = 0.5 as a mechanism of reinjection. 

The numerical results for maps described by Eqs. (HH, 
dni), dUl) and by Eqs. (US]), (|l5l)-(II71) are shown in Fig.[T] 
and Fig. [2l respectively. We calculate the power spectral 
density directly, according to the definition, as the nor- 
malized squared modulus of the Fourier transform of the 
signal. 



The numerical results for the map described by 
Eqs. dUl), dUl), daO]) are shown in Fig. |1 From Fig. [3i 
we can see that this map leads to the intermittent behav- 
ior. Due to larger average (z) and larger exponent rj the 
durations of laminar phases are shorter than in Figs. [T^, 
[2|a,. The PDF of the variable y, shown in Fig. [3Jd, has a 
power-law form with the exponent —3 for larger values 
of y, whereas for small values of y the PDF decreases ex- 
ponentially. The PSD of the time series {t/nji shown in 
Fig. [3t, has 1// behavior for a wide range of frequencies. 
The 1// interval in the PSD is 10"^ < / < lO'^. 

As the numerical examples show, both maps ([7|) and 
(jlOP for some values of the parameters can yield time 
series with 1// PSD in a wide range of frequencies. In 
addition, the PDF of the deviation from the invariant 
subspace y for these values of parameters has a power-law 
part with the exponent —3, in contrast to on-off intermit- 
tency where the exponent in the PDF is close to —1 and 
l/-y/7 PSD. The explanation of the observed behavior of 
PDF and PSD will be provided in the next Section. 



S{I) 



N 



(18) 



where the angle brackets (•) denote averaging over re- 
alizations. We used the time series {j/n} of the length 
N — 10^ and averaged over 100 realizations with ran- 
domly chosen initial value j/o- 

From Fig. [1^ and Fig. [2^ we can see that these maps 
indeed lead to intermittent behavior, where the laminar 
phases are changed by bursts of activity corresponding 
to the large deviations of the variable y from the aver- 
age value. The laminar phases of the first map appear 
smoother than laminar phases of the second. The PDF of 
the variable y, shown in Fig. [Hd and Fig. [2}d, has in both 
cases a power-law form with the exponent —3 for larger 
values of y, whereas for small values of y the PDF de- 
creases exponentially. The PSD of the time series {?/«}, 
shown in Fig.[TJ; and Fig. [It, has 1// behavior for a wide 
range of frequencies. The 1// interval in the PSD in 
Figs.[TJ:,[l: is 10"^ < / < lO""*. 

For the map (ITUl) we consider the case with 77 = 3. To 
avoid reaching of the limiting value s — 1/ [f] — 1) = 0.5 
we modify the map ^ by introducing the reflection from 
the boundary s„ = 0.5: 



Sn+l = 0.5- \S„ + Z„ 



0.51 



(19) 



Then the map (|10|) for the transformed variable y„ 
exp^(sn) takes the form 



yn+i 



1 

Vn 



zzji 



(20) 



For the variable a;„ we again use the tent map (|13|) and 
calculate z„ according to Eq. (fT4|) with the the average 
(z) = 9 X 10~^ and the variance((z — (z))^) = 1. 



IV. APPROXIMATION OF DISCRETE MAPS BY 
STOCHASTIC DIFFERENTIAL EQUATIONS 

To obtain analytical expressions for the PDF and PSD 
of the deviation y, we approximate the maps ([7]) and (jlOp 
by a SDE. To obtain the SDE corresponding to the map 
([7]) we proceed as follows: we replace the variable z„ by a 
random Gaussian variable having the same average and 
variance as z„ and interpret Eq. ([7]) as Euler-Marujama 
approximation of a SDE. In this way we get the following 
SDE: 



dy = a r/ 



V T] — I / y 



ymin \ \ 2»)-l 



y 



y^'^-'^dt+ay'^dW . 

(21) 
Here W^ is a standard Wiener process (the Brownian mo- 
tion) and the parameters a, 2/min, and v are given by the 
equations 



a=^{{z-{z)r), 
2(z) 



1^ = 2?] — 



(ry-l)((z-(z))2) 

27 



((Z-(Z))2) 



(22) 
(23) 
(24) 



SDE approximating the map (jlOp can be obtained in 
the following way: we approximate a random walk de- 
scribed by Eq. ([S]) by a Brownian motion with constant 
drift ds = adt + adW, where a is given by Eq. (l22t and 
a — (z). After transformation of the variable s to the 
variable y — exp (s) we get a particular case of Eq. (PT|) 
with V = rj, the other parameters a and ymin are given 
by Eqs. ^, ^. Thus, both maps © and ^ cor- 
respond to the same SDE ^. The SDE dH]) has the 
same form as that considered in Refs. 12. It is possible 
to obtain the non-linear SDE of the form 1^1^ starting 
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FIG. 1. (Color online) (a) Time series of map described by Eqs. (fT2)) . (fT3|) . (fTi)) . (b) PDF of the variable y. The dashed (green) 
line is the analytical expression (|29l) . (c) Power spectral density S{f) of time series. The dashed (green) line shows the slope 
1//. Parameters used are (z) = 5 x 10~^ , {{z ~ (z))"^) = 1. 



10 



(a) 



10"' 



ynio"^ 



10 



10 





k) 


WvAv- 



4 6 

n/lO** 



10 




FIG. 2. (Color onhne) (a) Time series of map described by Eqs. (fT3|l . (fT5|) - (fT7l) . (b) PDF of the variable y. The dashed (green) 
line is the analytical expression (|29|) . (c) Power spectral density S{f) of time series. The dashed (green) line shows the slope 
1//. Parameters used are the same as in Fig. [T] 



from the agent-based herding model.^^ In Ref. 28 modi- 
fications of these equations by introducing additional pa- 
rameters are presented. These equations may generate 
signals with q-exponential and q-Gaussian distribution 
of the nonextensive statistical mechanics. 

The approximation of the map ([7]) by the SDE ([2T1) is 
valid when the value of y is sufficiently small. The max- 
imum value of y can be determined from the condition 
that the second term in Eq. ([7]) should be much smaller 
than the first. We can estimate this condition as 



yinax 



V((z-(z))2)«y„ 



giving 



Vn 



< 



{{z-{z)r)-^^^ 



(25) 



(26) 



For the map ([TU]) the approximation by the SDE (PT|) 
is valid as long as the variable s„ is far from 1/(77 — 
1) where the g-exponential function exp (s„) becomes 
infinite. Assuming that the presence of the limiting value 
1/(77—1) does not infiuence the random walk ([9]) when the 
distance to this limiting value is larger than the standard 
deviation of Zn (that is, s„ < 1/(77 ~ 1) ~ "■) ^^ can 
estimate the maximum value of y as j/max S 6xp^ (1/(77 ~ 
I) — a). This estimation coincides with Eq. (P5|) . 

Using Eqs. (P5|) and (pS)) we get the expression for the 



ratio 77max/ymii 

yinax 
yniin 



< 



(77-l)v/((^-(z))2) 



2(z) 



(27) 



As it was shown in Refs. |12|, the SDE (|2T|) generates 
signals with power-law PSD in a wide range of frequencies 
when the variable y can vary in a wide region, j/max 3> 
T/min- The condition ymax/ymin > 1 is obeyed when 



((z-(z))2)»(z)^ 

that is, the standard deviation of the variable z„ 
be much larger than the average. 

The SDE JUl) leads to the steady state PDF 



(77- Dv"^^ 
^o(y) = ^V^T^exp 



r 



m 



y. 



Vn 



y 



v-i 



(28) 



should 



(29) 



Thus, the parameter v gives the exponent of the power- 
law part of the PDF and the parameter j/min gives the 
position of the exponential cut-off at small values of y. 
From Eq. (j23[) it follows that ymin grows with the growing 
average (z). As can be seen in Figs. [TJd, [5]d, [3]d, there is 
a good agreement of the numerically obtained PDF with 
the analytical expression (P5)l . Similarly as in the case of 
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FIG. 3. (Color online) (a) Time series of map described by Eqs. ITS]), ^T^, ((20l). (b) PDF of the variable y. The dashed (green) 
line is the analytical expression (|29l) . (c) Power spectral density S{f) of time series. The dashed (green) line shows the slope 
1//. Parameters used are (z) = 9 x 10""', {{z — (z))'^) = 1. 



on-ofF intermittency we obtain PDF of the deviation from 
the invariant subspace having power-law form, however, 
the power exponent v can assume values significantly dif- 
ferent from 1. 

Numerical analysis, performed in Ref.llSl indicates that 
the stochastic variable y, described by a SDE similar 
to (|2ip exhibits intermittent behavior: there are peaks, 
bursts or extreme events, corresponding to the large devi- 
ations of the variable from the appropriate average value, 
separated by laminar phases with a wide range distribu- 
tion of the laminar durations. The exponent —3/2 in the 
PDF of the interburst durations has been numerically 
obtained. 

In Refs. il^ it was shown that SDE (PTt generates sig- 
nals with PSD having the form S{f) ^ f^^ in a wide 
range of frequencies with the exponent 



/?=! 



I/-3 

2(^-1) 



(30) 



The connection of the PSD of the signal generated by 
SDE pT|) with the behavior of the eigenvalues of the 
corresponding Fokker-Planck equation was analyzed in 
Ref. Il4l An additional argument based on scaling prop- 
erties showing that PSD of the signal generated by 
SDE (|2T|1 has the power-law behavior in some range of 
frequencies we present in Appendix \^ For the parame- 
ters used in Figs. HJ HJ EJ Eq. §UI gives (3 = 1. Numer- 
ically obtained PSD shown in Figs. [TJ;, [5J:, [5J; confirms 
this prediction. Thus, as long as the approximation of 
the maps d?]) or (dUD by the SDE dU) is vahd, the PSD 
of the time series {yn} exhibits a power-law behavior, 
including 1// noise. 

The range of frequencies where PSD has power-law 
behavior is limited by the minimum and maximum values 
yinin and ymax- The limiting frequencies are estimated in 
Ref. [11 and also in Appendix[B] Using Eqs. ([22]), ((221) and 
(P^ . we can write the range of frequencies (jB4p where 
the PSD has the power-law form as 



Vn 



2(,,-l) 



If J/max/ymin ^ li this frequency range can span many 
orders of magnitude. However, this estimation of the 
frequency range is too broad and the numerical solution 
of Eq. (121]) gives much narrower range. Nevertheless, 
Eq. (|3T|) correctly reflects the following properties of the 
frequency region where PSD has 1//^ dependence: the 
width of this frequency region increases with increase of 
the ratio between minimum and maximum values, y„ 
and j/max, and with increase of the difference rj — Ir 



14 



< 27r/ < 1 



(31) 



V. CONCLUSIONS 

We demonstrate that the nonlinear maps having in- 
variant subspace and the expansion in the powers of the 
deviation from the invariant subspace having the form of 
Eq. (|4]) can generate signals with 1// noise. In contrast 
to known mechanism of 1// noise involving Pomeau- 
Manneville type maps, the parameter z„ in the map 
Eq. (HI is not static. Another difference is that the ex- 
ponent /? in the PSD, as Eq. ([30)) shows, depends on two 
parameters rj and v, thus 1//^ noise can be obtained for 
various values of the exponent /3. 

The width of the frequency region where the PSD 
has f~^ behavior is limited by the average value of the 
variable z„: this width increases as (z) approaches the 
threshold value (z) = 0. In addition, the width of the 
power-law region in the PSD increases with increasing 
the difference rj — 1. 



Appendix A: Nonlinear stochastic differential equation 
generating signals with 1//'* noise 

Pure 1/ f^ PSD is physically impossible because the 
total power would be infinity. Therefore we will consider 
signals with PSD having 1/ f^ behavior only in some wide 
intermediate region of frequencies, /min <C / ^ /max, 
whereas for small frequencies / <^ /min PSD is bounded. 
We can obtain nonlinear SDE generating signals exhibit- 
ing 1// noise using the following considerations. Wiener- 
Khintchine theorem relates PSD S{f) to the autocorre- 
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lation function C{t) 



S{f) cos{2Tr ft)dt. 



(Al) 



If S {f ) ^ f '^ina wide region of frequencies, then for the 
frequencies in this region the PSD has a scaling property 



S{af) ^ a-^Sif) 



(A2) 



when the influence of the limiting frequencies /min an 
/max is neglected. From the Wiener-Khintchine theorem 
(|A1[) it follows that the autocorrelation function has the 
scaling property 



C{at) - a'^~^C{t) 



(A3) 



in the time range l//max <^t <^ l//min- The autocorre- 
lation function can be written asi^i^i^S 

C{t) ^ J dy J dy' yy'Poiy)Py{y',t\y.,0) , (A4) 

where Po(y) is the steady-state PDF and Py{y' ,t\y,Q) 
is the transition probability (the conditional probability 
that at time t the signal has value y' with the condition 
that at time t = the signal had the value y). The 
transition probability can be obtained from the solution 
of the Fokker-Planck equation with the initial condition 
Py{y', t\y, 0) — 5[y' — y). The required property (|A3|) can 
be obtained when the steady-state PDF has the power- 
law form 



Po{y)^y-'' 



(A5) 



and the transition probability has the scaling property 

Py{ay', t\ay, 0) = a'^ Py{y' , a'^"-^'^, 0) , (A6) 

that is, change of the magnitude of the stochastic vari- 
able y is equivalent to the change of time scale. In this 
case from Eq. (|A4[) it follows that the autocorrelation 
function has the required property (jA3[) with /3 given by 
Eq. ([50]) . In order to avoid the divergence of steady state 
PDF (jA5[) the diffusion of stochastic variable y should 
be restricted at least from the side of small values and, 
therefore, Eq. (jASp holds only in some region of the 
variable y, ymin <C y ^ ymax- When the diffusion of 
stochastic variable y is restricted, Eq. (|A6|) also cannot 
be exact. However, if the influence of the limiting val- 
ues j/min and j/max Can be neglected for time t in some 
region tmin <C t ^ imax, we can expect that Eq. (|A3[) 
approximately holds for this time region. 

To get the required scaling (|A6p of the transition 
probability, the SDE should contain only powers of the 
stochastic variable y and the coefficient in the noise term 
should be proportional to y^ . The drift term then is fixed 
by the requirement (|A5[) for the steady-state PDF. Thus 
we consider SDE 



dy = a^[v- -^v 



f'^-^dt + ay'^dW . 



(A7) 



In order to obtain a stationary process and avoid the 
divergence of steady state PDF the diffusion of stochastic 
variable y should be restricted or equation (jA7|) should 
be modified. The simplest choice of the restriction is the 
reflective boundary conditions at y = j/min and y = ymax- 
Exponentially restricted diffusion with the steady state 
PDF 



Pa{y) ijexp 

V 



Vn 



V 



(A8) 



is generated by the SDE 



dy 



1 



2 V y- 



tfn 



y^'^-^dt+ay'^dW 

(A9) 

obtained from Eq. ()A7|) by introducing the additional 
terms. 



Appendix B: Estimation of the range of the frequencies 
where PSD has the power-law behavior 

The presence of the restrictions at y = j/min and 
y = ymax makes the scaling (jA6|) not exact and this lim- 
its the power-law part of the PSD to a finite range of 
frequencies /min <C / ^ /max- Let us estimate the limit- 
ing frequencies. Taking into account the limiting values 
j/min and ymax, Eq. (JA6| for the transition probability 
corresponding to SDE (jA7l) becomes 

Py(ay',i|ay, 0;aymin,a2/max) = 

a~^Py{y', a^^''^~^h\y, 0; ymi„, ymax) • (Bl) 

The steady-state distribution Po(y;2/min,ymax) has the 
scaling property 

Pa{ay; aymin, aymax) = a^^Poiy, ymin, ymax) • (B2) 
Inserting Eqs. (IBTJI and ((B2|) into Eq. (|X4|) we obtain 

C{t; aymin, aymax) = a^C{a^^'''^h, ymin, ymax) ■ (B3) 

This equation means that time t in the autocorrelation 

function should enter only in combinations with the lim- 

1 1 

iting values, ymin^^'""^' and ymax*^*''"^' ■ We can expect 

that the influence of the limiting values can be neglected 

and Eq. (jA6p holds when the first combination is small 

and the second large, that is when time t is in the interval 

fT-^y^^i-") < t < (T-^y^^i-"). Then, using Eq. ^^ the 

frequency range where the PSD has 1/ f^ behavior can 

be estimated as 



a2y^([; 1' « 2^/ « a'y 
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